Zygotic Splicing Activation of the Transcriptome is a Crucial Aspect of Maternal-to-Zygotic Transition and Required for the Conversion from Totipotency to Pluripotency.

During maternal-to-zygotic transition (MZT) in the embryo, mRNA undergoes complex post-transcriptional regulatory processes. However, it is unclear whether and how alternative splicing plays a functional role in MZT. By analyzing transcriptome changes in mouse and human early embryos, dynamic changes in alternative splicing during MZT are observed and a previously unnoticed process of zygotic splicing activation (ZSA) following embryonic transcriptional activation is described. As the underlying mechanism of RNA splicing, splicing factors undergo dramatic maternal-to-zygotic conversion. This conversion relies on the key maternal factors BTG4 and PABPN1L and is zygotic-transcription-dependent. CDK11-dependent phosphorylation of the key splicing factor, SF3B1, and its aggregation with SRSF2 in the subnuclear domains of 2-cell embryos are prerequisites for ZSA. Isoforms generated by erroneous splicing, such as full-length Dppa4, hinder normal embryonic development. Moreover, alternative splicing regulates the conversion of early embryonic blastomeres from totipotency to pluripotency, thereby affecting embryonic lineage differentiation. ZSA is an essential post-transcriptional process of MZT and has physiological significance in generating new life. In addition to transcriptional activation, appropriate expression of transcript isoforms is also necessary for preimplantation embryonic development.

Different microbial communities in paddy soils under organic and nonorganic farming.

Organic agriculture is a farming method that provides healthy food and is friendly to the environment, and it is developing rapidly worldwide. This study compared microbial communities in organic farming (Or) paddy fields to those in nonorganic farming (Nr) paddy fields based on 16S rDNA sequencing and analysis. The predominant microorganisms in both soils were Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria, and Nitrospirota. The alpha diversity of the paddy soil microbial communities was not different between the nonorganic and organic farming systems. The beta diversity of nonmetric multidimensional scaling (NMDS) revealed that the two groups were significantly separated. Distance-based redundancy analysis (db-RDA) suggested that soil pH and electrical conductivity (EC) had a positive relationship with the microbes in organic paddy soils. There were 23 amplicon sequence variants (ASVs) that showed differential abundance. Among them, g_B1-7BS (Proteobacteria), s_Sulfuricaulis limicola (Proteobacteria), g_GAL15 (p_GAL15), c_Thermodesulfovibrionia (Nitrospirota), two of f_Anaerolineaceae (Chloroflexi), and two of g_S085 (Chloroflexi) showed that they were more abundant in organic soils, whereas g_11-24 (Acidobacteriota), g__Subgroup_7 (Acidobacteriota), and g_Bacillus (Firmicutes) showed differential abundance in nonorganic paddy soils. Functional prediction of microbial communities in paddy soils showed that functions related to carbohydrate metabolism could be the major metabolic activities. Our work indicates that organic farming differs from nonorganic farming in terms of microbial composition in paddy soils and provides specific microbes that might be helpful for understanding soil fertility.

Genome-wide map of R-loops reveals its interplay with transcription and genome integrity during germ cell meiosis.

INTRODUCTION: The R-loop is a naturally formed three-strand nucleic acid structure that recently has been reported to participate in multiple biological processes and helped answer some previously unexplained scientific questions. Meiosis process involves multiple chromatin-related events such as DNA double-stranded breaks (DSB) formation, repairing and transcriptional dynamics. OBJECTIVES: Explore the regulatory roles and physiological functions of R-loops in the mammalian meiosis process. METHODS: In our study, using genome-wide S9.6 CUT & Tag seq, we first mapped the genomic distribution and dynamic changes of R-loop during the meiotic process in mice, from spermatogonia to secondary spermatocytes. And we further explore the role of R-loop in physiological conditions by constructing conditional knockout mice of Rnaseh1, which deleted the R-loop endonuclease before meiosis entry. RESULTS: R-loop predominantly distributes at promoter-related regions and varies across different meiotic stages. By joint analysis with the corresponding transcriptome, we found that the R-loop was closely related to transcription during the meiotic process. The high frequency of promoter-related R-loop in meiotic cells is usually accompanied by high transcription activity, and we further verified this in the leptotene/zygotene to the pachytene transition process. Moreover, the lack of RNase H1 caused sterility in male mice with R-loop accumulation and abnormal DSB repair in spermatocytes. Further analysis showed that abnormal R-loop accumulation in the leptotene/zygotene stages influenced transcriptional regulation in the pachytene stage. CONCLUSION: The mutual regulation of the R-loop and transcription plays an essential role in spermatogenesis. And R-loop is also important for the normal repair process of DSB during meiosis.

PABPN1 functions as a hub in the assembly of nuclear poly(A) domains that are essential for mouse oocyte development.

Growing oocytes store a large amount of maternal mRNA to support the subsequent "maternal-zygotic transition" process. At present, it is not clear how the growing oocytes store and process the newly transcribed mRNA under physiological conditions. In this study, we report non-membrane-bound compartments, nuclear poly(A) domains (NPADs), as the hub for newly transcribed mRNA, in developing mouse oocytes. The RNA binding protein PABPN1 promotes the formation of NPAD through its N-terminal disordered domain and RNA-recognized motif by means of liquid phase separation. Pabpn1-null growing oocytes cannot form NPAD normally in vivo and have defects in stability of oocyte growing-related transcripts and formation of long 3' untranslated region isoform transcripts. Ultimately, Pabpn1fl/fl;Gdf9-Cre mice are completely sterile with primary ovarian insufficiency. These results demonstrate that NPAD formed by the phase separation properties of PABPN1-mRNA are the hub of the newly transcribed mRNA and essential for the development of oocytes and female reproduction.

The RNA ligase RNA terminal phosphate cyclase B regulates mRNA alternative splicing and is required for mouse oocyte development and maintenance.

Recent large-scale mRNA sequencing has shown that introns are retained in 5-10% of mRNA, and these events are named intron retention (IR). IR has been recognized as a key mechanism in the regulation of gene expression. However, the role of this mechanism in female reproduction in mammals remains unclear. RNA terminal phosphate cyclase B (RTCB) is a RNA ligase; we found that RTCB conditional knockout mice have premature ovarian failure and that RTCB plays a crucial role in follicular development. RTCB regulated the splicing of transcripts related to DNA methylation and DNA damage repair. In addition, it regulated the resumption of oocyte meiosis by affecting CDK1 activation. Moreover, the loss of RTCB suppressed zygotic genome activation (ZGA) and decreased translation at the global level. In addition, Rtcb deletion resulted in the accumulation of maternal mRNAs containing unspliced introns and in a decline in the overall level of transcripts. As a result, the Rtcb-/- females were sterile. Our study highlights the important role of RTCB-regulated noncanonical alternative splicing in female reproduction.

NAT10-mediated N4-acetylcytidine modification is required for meiosis entry and progression in male germ cells.

Post-transcriptional RNA modifications critically regulate various biological processes. N4-acetylcytidine (ac4C) is an epi-transcriptome, which is highly conserved in all species. However, the in vivo physiological functions and regulatory mechanisms of ac4C remain poorly understood, particularly in mammals. In this study, we demonstrate that the only known ac4C writer, N-acetyltransferase 10 (NAT10), plays an essential role in male reproduction. We identified the occurrence of ac4C in the mRNAs of mouse tissues and showed that ac4C undergoes dynamic changes during spermatogenesis. Germ cell-specific ablation of Nat10 severely inhibits meiotic entry and leads to defects in homologous chromosome synapsis, meiotic recombination and repair of DNA double-strand breaks during meiosis. Transcriptomic profiling revealed dysregulation of functional genes in meiotic prophase I after Nat10 deletion. These findings highlight the crucial physiological functions of ac4C modifications in male spermatogenesis and expand our understanding of its role in the regulation of specific physiological processes in vivo.

Dynamic mRNA degradome analyses indicate a role of histone H3K4 trimethylation in association with meiosis-coupled mRNA decay in oocyte aging.

A decrease in oocyte developmental potential is a major obstacle for successful pregnancy in women of advanced age. However, the age-related epigenetic modifications associated with dynamic transcriptome changes, particularly meiotic maturation-coupled mRNA clearance, have not been adequately characterized in human oocytes. This study demonstrates a decreased storage of transcripts encoding key factors regulating the maternal mRNA degradome in fully grown oocytes of women of advanced age. A similar defect in meiotic maturation-triggered mRNA clearance is also detected in aged mouse oocytes. Mechanistically, the epigenetic and cytoplasmic aspects of oocyte maturation are synchronized in both the normal development and aging processes. The level of histone H3K4 trimethylation (H3K4me3) is high in fully grown mouse and human oocytes derived from young females but decreased during aging due to the decreased expression of epigenetic factors responsible for H3K4me3 accumulation. Oocyte-specific knockout of the gene encoding CxxC-finger protein 1 (CXXC1), a DNA-binding subunit of SETD1 methyltransferase, causes ooplasm changes associated with accelerated aging and impaired maternal mRNA translation and degradation. These results suggest that a network of CXXC1-maintained H3K4me3, in association with mRNA decay competence, sets a timer for oocyte deterioration and plays a role in oocyte aging in both mouse and human oocytes.

USP16-mediated histone H2A lysine-119 deubiquitination during oocyte maturation is a prerequisite for zygotic genome activation.

Maternal-to-zygotic transition (MZT) is the first and key step in the control of animal development and intimately related to changes in chromatin structure and histone modifications. H2AK119ub1, an important epigenetic modification in regulating chromatin configuration and function, is primarily catalyzed by PRC1 and contributes to resistance to transcriptional reprogramming in mouse embryos. In this study, the genome-wide dynamic distribution of H2AK119ub1 during MZT in mice was investigated using chromosome immunoprecipitation and sequencing. The results indicated that H2AK119ub1 accumulated in fully grown oocytes and was enriched at the TSSs of maternal genes, but was promptly declined after meiotic resumption at genome-wide including the TSSs of early zygotic genes, by a previously unidentified mechanism. Genetic evidences indicated that ubiquitin-specific peptidase 16 (USP16) is the major deubiquitinase for H2AK119ub1 in mouse oocytes. Conditional knockout of Usp16 in oocytes did not impair their survival, growth, or meiotic maturation. However, oocytes lacking USP16 have defects when undergoing zygotic genome activation or gaining developmental competence after fertilization, potentially associated with high levels of maternal H2AK119ub1 deposition on the zygotic genomes. Taken together, H2AK119ub1 level is declined during oocyte maturation by an USP16-dependent mechanism, which ensures zygotic genome reprogramming and transcriptional activation of essential early zygotic genes.

Nuclear poly(A) binding protein 1 (PABPN1) mediates zygotic genome activation-dependent maternal mRNA clearance during mouse early embryonic development.

An embryo starts its life with maternal mRNA clearance, which is crucial for embryonic development. The elimination of maternal transcripts occurs by the joint action of two pathways: the maternally encoded mRNA decay pathway (M-decay) and the zygotic genome activation (ZGA)-dependent pathway (Z-decay). However, zygotic factors triggering maternal mRNA decay in early mammalian embryos remain largely unknown. In this study, we identified the zygotically encoded nuclear poly(A) binding protein 1 (PABPN1) as a factor required for maternal mRNA turnover, with a previously undescribed cytoplasmic function. Cytoplasmic PABPN1 docks on 3'-uridylated transcripts, downstream of terminal uridylyl transferases TUT4 and TUT7, and recruits 3'-5' exoribonuclease DIS3L2 to its targets, facilitating maternal mRNA decay. Pabpn1-knockout in mice resulted in preimplantation stage mortality due to early developmental arrest at the morula stage. Maternal mRNAs to be eliminated via the Z-decay pathway failed to be removed from Pabpn1-depleted embryos. Furthermore, PABPN1-mediated Z-decay is essential for major ZGA and regulates the expression of cell fate-determining factors in mouse preimplantation embryos. This study revealed an unforeseen cytoplasmic function of PABPN1 coupled with early embryonic development, characterized the presence of a zygotic destabilizer of maternal mRNA, and elucidated the Z-decay process mechanisms, which potentially contribute to human fertility.

CNOT6/6L-mediated mRNA degradation in ovarian granulosa cells is a key mechanism of gonadotropin-triggered follicle development.

CCR4-NOT deadenylase is a major regulator of mRNA turnover. It contains two heterogeneous catalytic subunits CNOT7/8 and CNOT6/6L in vertebrates. The physiological function of each catalytic subunit is unclear due to the gene redundancy. In this study, Cnot6/6l double knockout mice are generated. Cnot6l-/- female mice are infertile, with poor ovarian responses to gonadotropins. Follicle-stimulating hormone (FSH) stimulates the transcription and translation of Cnot6 and Cnot6l in ovarian granulosa cells. CNOT6/6L function as key effectors of FSH in granulosa cells and trigger the clearance of specific transcripts in granulosa cells during preantral to antral follicle transition. These results demonstrate that FSH modulates granulosa cell function by stimulating selective translational activation and degradation of existing mRNAs, in addition to inducing de novo gene transcription. Meanwhile, this study provides in vivo evidence that CNOT6/6L-mediated mRNA deadenylation is dispensable in most somatic cell types, but is essential for female reproductive endocrine regulation.

The CNOT4 Subunit of the CCR4-NOT Complex is Involved in mRNA Degradation, Efficient DNA Damage Repair, and XY Chromosome Crossover during Male Germ Cell Meiosis.

The CCR4-NOT complex is a major mRNA deadenylase in eukaryotes, comprising the catalytic subunits CNOT6/6L and CNOT7/8, as well as CNOT4, a regulatory subunit with previously undetermined functions. These subunits have been hypothesized to play synergistic biochemical functions during development. Cnot7 knockout male mice have been reported to be infertile. In this study, viable Cnot6/6l double knockout mice are constructed, and the males are fertile. These results indicate that CNOT7 has CNOT6/6L-independent functions in vivo. It is also demonstrated that CNOT4 is required for post-implantation embryo development and meiosis progression during spermatogenesis. Conditional knockout of Cnot4 in male germ cells leads to defective DNA damage repair and homologous crossover between X and Y chromosomes. CNOT4 functions as a previously unrecognized mRNA adaptor of CCR4-NOT by targeting mRNAs to CNOT7 for deadenylation of poly(A) tails, thereby mediating the degradation of a subset of transcripts from the zygotene to pachytene stage. The mRNA removal promoted by the CNOT4-regulated CCR4-NOT complex during the zygotene-to-pachytene transition is crucial for the appropriate expression of genes involved in the subsequent events of spermatogenesis, normal DNA double-strand break repair during meiosis, efficient crossover between X and Y chromosomes, and ultimately, male fertility.

Role of CxxC-finger protein 1 in establishing mouse oocyte epigenetic landscapes.

During oogenesis, oocytes gain competence and subsequently undergo meiotic maturation and prepare for embryonic development; trimethylated histone H3 on lysine-4 (H3K4me3) mediates a wide range of nuclear events during these processes. Oocyte-specific knockout of CxxC-finger protein 1 (CXXC1, also known as CFP1) impairs H3K4me3 accumulation and causes changes in chromatin configurations. This study investigated the changes in genomic H3K4me3 landscapes in oocytes with Cxxc1 knockout and the effects on other epigenetic factors such as the DNA methylation, H3K27me3, H2AK119ub1 and H3K36me3. H3K4me3 is overall decreased after knocking out Cxxc1, including both the promoter region and the gene body. CXXC1 and MLL2, which is another histone H3 methyltransferase, have nonoverlapping roles in mediating H3K4 trimethylation during oogenesis. Cxxc1 deletion caused a decrease in DNA methylation levels and affected H3K27me3 and H2AK119ub1 distributions, particularly at regions with high DNA methylation levels. The changes in epigenetic networks implicated by Cxxc1 deletion were correlated with the transcriptional changes in genes in the corresponding genomic regions. This study elucidates the epigenetic changes underlying the phenotypes and molecular defects in oocytes with deleted Cxxc1 and highlights the role of CXXC1 in orchestrating multiple factors that are involved in establishing the appropriate epigenetic states of maternal genome.

Genomewide decoupling of H2AK119ub1 and H3K27me3 in early mouse development.

Polycomb group (PcG) proteins are crucial chromatin regulators during development. H2AK119ub1 (H2Aub) and H3K27me3 are catalyzed by Polycomb-repressive complex 1 and 2 (PRC1/2) respectively, and they largely overlap in the genome due to mutual recruitment of the two complexes. However, it is unclear whether PRC1/H2Aub and PRC2/H3K27me3 can also function independently. By developing an ultra-sensitive carrier-DNA-assisted chromatin immunoprecipitation sequencing method termed CATCH-Seq, we generated allelic H2Aub profiles in mouse gametes and early embryos. Our results revealed an unexpected genomewide decoupling of H2Aub and H3K27me3 in mouse preimplantation embryos, where H2Aub but not H3K27me3 was enriched at PcG targets while only H3K27me3 was deposited in the broad distal domains associated with DNA methylation-independent non-canonical imprinting. These observations suggest that H2Aub represses future bivalent genes during early embryogenesis without H3K27me3, but it is not required for the maintenance of non-canonical imprinting, which is mediated by maternal H3K27me3. Thus, our study reveals the distinct depositions and independent functions of H2Aub and H3K27me3 during early mammalian development.

Dynamics and clinical relevance of maternal mRNA clearance during the oocyte-to-embryo transition in humans.

Maternal mRNA clearance is an essential process that occurs during maternal-to-zygotic transition (MZT). However, the dynamics, functional importance, and pathological relevance of maternal mRNA decay in human preimplantation embryos have not yet been analyzed. Here we report the zygotic genome activation (ZGA)-dependent and -independent maternal mRNA clearance processes during human MZT and demonstrate that subgroups of human maternal transcripts are sequentially removed by maternal (M)- and zygotic (Z)-decay pathways before and after ZGA. Key factors regulating M-decay and Z-decay pathways in mouse have similar expression pattern during human MZT, suggesting that YAP1-TEAD4 transcription activators, TUT4/7-mediated mRNA 3'-oligouridylation, and BTG4/CCR4-NOT-induced mRNA deadenylation may also be involved in the regulation of human maternal mRNA stability. Decreased expression of these factors and abnormal accumulation of maternal transcripts are observed in the development-arrested embryos of patients who seek assisted reproduction. Defects of M-decay and Z-decay are detected with high incidence in embryos that are arrested at the zygote and 8-cell stages, respectively. In addition, M-decay is not found to be affected by maternal TUBB8 mutations, although these mutations cause meiotic cell division defects and zygotic arrest, which indicates that mRNA decay is regulated independent of meiotic spindle assembly. Considering the correlations between maternal mRNA decay defects and early developmental arrest of in vitro fertilized human embryos, M-decay and Z-decay pathway activities may contribute to the developmental potential of human preimplantation embryos.

PABPN1L mediates cytoplasmic mRNA decay as a placeholder during the maternal-to-zygotic transition.

Maternal mRNA degradation is a critical event of the maternal-to-zygotic transition (MZT) that determines the developmental potential of early embryos. Nuclear Poly(A)-binding proteins (PABPNs) are extensively involved in mRNA post-transcriptional regulation, but their function in the MZT has not been investigated. In this study, we find that the maternally expressed PABPN1-like (PABPN1L), rather than its ubiquitously expressed homolog PABPN1, acts as an mRNA-binding adapter of the mammalian MZT licensing factor BTG4, which mediates maternal mRNA clearance. Female Pabpn1l null mice produce morphologically normal oocytes but are infertile owing to early developmental arrest of the resultant embryos at the 1- to 2-cell stage. Deletion of Pabpn1l impairs the deadenylation and degradation of a subset of BTG4-targeted maternal mRNAs during the MZT. In addition to recruiting BTG4 to the mRNA 3'-poly(A) tails, PABPN1L is also required for BTG4 protein accumulation in maturing oocytes by protecting BTG4 from SCF-ßTrCP1 E3 ubiquitin ligase-mediated polyubiquitination and degradation. This study highlights a noncanonical cytoplasmic function of nuclear poly(A)-binding protein in mRNA turnover, as well as its physiological importance during the MZT.

Genome-wide survey of the bHLH super gene family in Brassica napus.

BACKGROUND: The basic helix-loop-helix (bHLH) gene family is one of the largest transcription factor families in plants and is functionally characterized in diverse species. However, less is known about its functions in the economically important allopolyploid oil crop, Brassica napus. RESULTS: We identified 602 potential bHLHs in the B. napus genome (BnabHLHs) and categorized them into 35 subfamilies, including seven newly separated subfamilies, based on phylogeny, protein structure, and exon-intron organization analysis. The intron insertion patterns of this gene family were analyzed and a total of eight types were identified in the bHLH regions of BnabHLHs. Chromosome distribution and synteny analyses revealed that hybridization between Brassica rapa and Brassica oleracea was the main expansion mechanism for BnabHLHs. Expression analyses showed that BnabHLHs were widely in different plant tissues and formed seven main patterns, suggesting they may participate in various aspects of B. napus development. Furthermore, when roots were treated with five different hormones (IAA, auxin; GA3, gibberellin; 6-BA, cytokinin; ABA, abscisic acid and ACC, ethylene), the expression profiles of BnabHLHs changed significantly, with many showing increased expression. The induction of five candidate BnabHLHs was confirmed following the five hormone treatments via qRT-PCR. Up to 246 BnabHLHs from nine subfamilies were predicted to have potential roles relating to root development through the joint analysis of their expression profiles and homolog function. CONCLUSION: The 602 BnabHLHs identified from B. napus were classified into 35 subfamilies, and those members from the same subfamily generally had similar sequence motifs. Overall, we found that BnabHLHs may be widely involved in root development in B. napus. Moreover, this study provides important insights into the potential functions of the BnabHLHs super gene family and thus will be useful in future gene function research.

ZAR1 and ZAR2 are required for oocyte meiotic maturation by regulating the maternal transcriptome and mRNA translational activation.

Zar1 was one of the earliest mammalian maternal-effect genes to be identified. Embryos derived from Zar1-null female mice are blocked before zygotic genome activation; however, the underlying mechanism remains unclear. By knocking out Zar1 and its homolog Zar2 in mice, we revealed a novel function of these genes in oocyte meiotic maturation. Zar1/2-deleted oocytes displayed delayed meiotic resumption and polar body-1 emission and a higher incidence of abnormal meiotic spindle formation and chromosome aneuploidy. The grown oocytes of Zar1/2-null mice contained decreased levels of many maternal mRNAs and displayed a reduced level of protein synthesis. Key maturation-associated changes failed to occur in the Zar1/2-null oocytes, including the translational activation of maternal mRNAs encoding the cell-cycle proteins cyclin B1 and WEE2, as well as maternal-to-zygotic transition (MZT) licensing factor BTG4. Consequently, maternal mRNA decay was impaired and MZT was abolished. ZAR1/2 bound mRNAs to regulate the translational activity of their 3'-UTRs and interacted with other oocyte proteins, including mRNA-stabilizing protein MSY2 and cytoplasmic lattice components. These results countered the traditional view that ZAR1 only functions after fertilization and highlight a previously unrecognized role of ZAR1/2 in regulating the maternal transcriptome and translational activation in maturing oocytes.

Evolutionary and Comparative Expression Analyses of TCP Transcription Factor Gene Family in Land Plants.

The plant-specific Teosinte-branched 1/Cycloidea/Proliferating (TCP) transcription factor genes are involved in plants' development, hormonal pathways, and stress response but their evolutionary history is uncertain. The genome-wide analysis performed here for 47 plant species revealed 535 TCP candidates in terrestrial plants and none in aquatic plants, and that TCP family genes originated early in the history of land plants. Phylogenetic analysis divided the candidate genes into Classes I and II, and Class II was further divided into CYCLOIDEA (CYC) and CINCINNATA (CIN) clades; CYC is more recent and originated from CIN in angiosperms. Protein architecture, intron pattern, and sequence characteristics were conserved in each class or clade supporting this classification. The two classes significantly expanded through whole-genome duplication during evolution. Expression analysis revealed the conserved expression of TCP genes from lower to higher plants. The expression patterns of Class I and CIN genes in different stages of the same tissue revealed their function in plant development and their opposite effects in the same biological process. Interaction network analysis showed that TCP proteins tend to form protein complexes, and their interaction networks were conserved during evolution. These results contribute to further functional studies on TCP family genes.

Global Analysis of WOX Transcription Factor Gene Family in Brassica napus Reveals Their Stress- and Hormone-Responsive Patterns.

The plant-specific WUSCHEL-related homeobox (WOX) transcription factor gene family is important for plant growth and development but little studied in oil crops. We identified and characterized 58 putative WOX genes in Brassica napus (BnWOXs), which were divided into three major clades and nine subclades based on the gene structure and conserved motifs. Collinearity analysis revealed that most BnWOXs were the products of allopolyploidization and segmental duplication events. Gene structure analysis indicated that introns/exons and protein motifs were conserved in each subclade and RNA sequencing revealed that BnWOXs had narrow expression profiles in major tissues and/or organs across different developmental stages. The expression pattern of each clade was highly conserved and similar to that of the sister and orthologous pairs from Brassica rapa and Brassica oleracea. Quantitative real-time polymerase chain reaction showed that members of the WOX4 subclade were induced in seedling roots by abiotic and hormone stresses, indicating their contribution to root development and abiotic stress responses. 463 proteins were predicted to interact with BnWOXs, including peptides regulating stem cell homeostasis in meristems. This study provides insights into the evolution and expression of the WOX gene family in B. napus and will be useful in future gene function research.